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In the present work we propose a theory for obtaining successively better approximations to 
the linear response functions of time-dependent density or current-density functional theory. The 
new technique is based on the variational approach to many-body perturbation theory (MBPT) 
as developed during the sixties and later expanded by us in the mid nineties. Due to this feature 
the resulting response functions obey a large number of conservation laws such as particle and 
momentum conservation and sum rules. The quality of the obtained results is governed by the 
physical processes built in through MBPT but also by the choice of variational expressions. We here 
present several conserving response functions of different sophistication to be used in the calculation 
of the optical response of solids and nano-scale systems. 

PACS numbers: 



I. INTRODUCTION 



Optical spectra constitute important tools for gain- 
ing information on the electronic structure of solids, 
molecules, and nano-systems. In many systems the 
particle-hole interaction leads to a strong excitonic dis- 
tortion of the optical spectrum - particularly in nano- 
scale objects. The theoretical description of such spectra 
is relatively sophisticated and very costly from a com- 
putational point of view. Some time ago, it was real- 
ized that these spectra are also within reach using time- 
dependent (TD) density-functional theory (DFT) - but 
with much less computational effort. From the Runge- 
Gross theorem-! of TDDFT we know how to construct the 
exact density response function of any electronic system 
in terms of an exchange-correlation kernel describing the 
particle-hole interactions. And from recent work by sev- 
eral researchers^i^ we have a rather good idea about the 
properties of this kernel if it is to reproduce the rather ac- 
curate results obtained from solving the Bethe-Salpeter 
equation of many-body perturbation theory (MBPT). 
The kernel, usually named /xc, has been calculated in 
the exchange-only approximation of TDDFT by several 
people in the past, see for instance Refs. 00. More re- 
cently, the kernel /xc has been calculated in the same 
approximation by Petersilka, Gossmann, and Gross^ for 
the helium atom, by Kurth and von Barth for the den- 
sity response of the homogeneous electron gas^ and by 
Kim and Gorling^ in the case of bulk silicon. In the cases 
of atomic helium and the homogeneous electron gas the 
resulting response function represented a substantial im- 
provement on that of the Random Phase Approximation 
(RPA). The excitation energies of helium were much im- 
proved and the total energies obtained from the response 
function were much superior to those obtained from the 
RPA response function in both helium and the homoge- 
neous electron gas. Unfortunately, this ab initio approach 
did not work very well in bulk silicon unless one rather 
arbitrarily introduces some kind of static screening of the 



particle-hole interaction. 

In actual fact, within TDDFT, no systematic and re- 
alistic route toward successively better approximations 
has, so far, been available. In the present work we have 
constructed such a scheme based on the variational ap- 
proach to many-body theory developed in Ref 11. In 
terms of the one-electron Green function of MBPT, these 
functionals give stationary expressions for the total ac- 
tion of the system at hand - or the total energy in the case 
of time independent problems. From a stationary action 
it is rather straight-forward to construct the time depen- 
dent density response function. Building the functionals 
from the ^-derivable theory of Baym and KadanoffjiSii^ 
always results in response functions which obey essential 
physical constraints like particle, momentum or energy 
conservation. 

The simple idea of the present work is to restrict the 
variational freedom of the functionals to the domain of 
Green functions which are non-interacting and given by 
a local one-electron potential - and vector-potential in 
case of current-DFT. According to the Runge-Gross the- 
orem this restriction immediately results in a density- 
functional theory the quality of which is determined by 
the sophistication which is build in to the choice of <I> 
derivable approximation for the action functional. Thus, 
to every conserving scheme within MBPT there is a cor- 
responding level of approximation within TDDFT. The 
latter is determined variationally and there is no longer 
a need for an ad hoc procedure to equate corresponding 
quantities between TDDFT and MBPT. A potentially in- 
teresting consequence of the theory proposed here, is that 
the often discussed linearized Sham-Schliiter equationi^ 
for the exchange-correlation potential is nothing but the 
stationary condition for the action functional. In the par- 
ticular version of the variational functionals developed 
in Ref. and named 4" derivable theories, also the 
screened Coulomb interaction becomes an independent 
variable at ones disposal. This leads to approximations 
within TDDFT which are potentially as accurate as those 
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of more elaborate schemes within MBPT but which are 
comparatively easier to implement - especially in nano- 
systems and complex solids. 

II. VARIATIONAL APPROACH TO TDDFT 

Let us consider a system of interacting fermions ex- 
posed to an external, possibly time-dependent field 
w{rt). The full many-body Hamiltonian reads 

H = f + U + W, (1) 

where 

is the kinetic operator, while 

(fird^r' i;^r)i;^{r')v{r,r')^P{r')^P{r), 

is the interaction operator {v{r,r') = l/\r — r'\). The 
coupling to the external field is given by 

= y d^r w{rt)n{r), 

where n{r) = ip^{r)ip{r) is the density operator. The 
Green function G obeys Dyson's equation 

G — Gfj + GhSG 

where Gh is the Hartree Green function and E is the 
exchange-correlation part of the electronic self-energy. 
Diagrammatic perturbation theory provides a tool for 
generating approximate self-energies and, in turn, ap- 
proximate Green functions. Except for physical intu- 
ition, the diagrammatic techniques rely solely on the 
validity of Wick's theorem, -^-ifi Thus, a typical contri- 
bution to the self-energy is represented by a diagram 
containing non-interacting propagators and interaction 
lines. However, any approximation which contains only 
a finite number of these diagrams violates many conser- 
vation laws. Conserving approximations require a proper 
choice of an infinite set of diagrams. The conserving 
approach by Baym~- was based on such choices. Also 
the variational scheme by Almbladh, von Barth and van 
Leeuwenii (ABL) was designed with the same objective 
in mind. The former approach is referred to as a $- 
derivable scheme because its central quantity is a univer- 
sal functional, called <I>, of the one-electron Green func- 
tion G and the bare Coulomb potential v. It is con- 
structed such that its functional derivative with respect 
to G gives the exchange-correlation part of the electronic 
self-energy S whereas the functional derivative with re- 
spect to the Coulomb interaction v essentially gives the 
reducible polarizability x of the system. 



(Here and in the following we use the short-hand nota- 
tion 1 = (rijfi), 2 = (r2,i2) and so on). Notice, how- 
ever, that there is no reference to an actual system in 
the functional. It acquires a meaning only when it is 
evaluated at a Green function of an actual system. In the 
approach of ABL, the central quantity is instead the func- 
tional 'J having the Green function G and the screened 
Coulomb interaction W as independent variables. It is 
constructed so as to give the self-energy when it is dif- 
ferentiated with respect to G and the irreducible polar- 
izability P when functionally differentiated with respect 
to W. Again, there is no reference to the actual sys- 
tem contained in the functional 4". By adding functional 
pieces to the 4> or the ^' functional respectively, pieces 
which do contain clear connections to the system under 
study (like, e.g., the externally applied potential w), one 
constructs functionals for the total energy - or the ac- 
tion in the case of time dependent problems - which, as 
functionals of G, have their stationary point at the Green 
function G which is the solution to Dyson's equation. In 
the case of the Vt-based functionals they are also station- 
ary when the screened interaction W obeys the so called 
reduced Bethe-Salpeter equation to be discussed later. 

The first variational functional of this kind was con- 
structed by Luttinger and Wardii (LW). It is a $ func- 
tional and it has the appearance 

iYlw[G] = $[G] - Tr {EG + ln(E - G^^)} - iUaiG]. 

(3) 

In Eq. ©, the functional Uh[G] = -^Tt{VhG} is 
the classical Hartree energy, Vnirt) ~ J d^r'v{Y, r')n{r't) 
is the Hartree potential and n{rt) the electron density. 
The symbol Tr (Trace) denotes a sum over labels of 
one-electron states plus an integration over time, or fre- 
quency for equilibrium problems in the ground state, a 
sum over discrete frequencies or an integral over imagi- 
nary times for elevated temperatures or an integral along 
the Keldysh contouriiSiiSiiS in the case of non-equilibrium 
problemsiSS It is straightforward to realize that Ilw is 
stationary when G obeys Dyson's equation with the self- 
energy of Eq. ||2Jl. At the stationary point, the Green 
function is fully conserving. 

At this point we would like to draw the readers atten- 
tion to a very interesting fact, the ramifications of which 
have yet to be discovered. The variational schemes are, 
by no means, unique. By adding to YlWi any functional 
F[D], where 

D[G] = G(Gh^ - E[G]) - 1, (4) 

obeying 

one obtains a new variational functional having the same 
stationary point and the same value at the stationary 
point. It might, however, be designed to give a second 
derivative which also vanishes at the stationary point - 
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something that would be of utmost practical value. Such 
possibilities could open up a whole new field of research. 
Choosing to add F[D] to the LW functional, where 



F[D] = Tr {-D + lii{D + 1)} 



(5) 



obviously has the desired properties, leads to the func- 
tional 

iYK[G] = $[G] - Tr {GGh^ - 1 + ln(-G"^)} - iUniG]. 

This functional was first written down by Klein, and 
could thus be called the Klein functional in order to 
distinguish it from the LW functional above. Unfortu- 
nately, this functional is less stable (large second deriva- 
tive) at the stationary point as compared to the LW 
functional Since the construction of response 
functions for TDDFT from the variational functionals 
involve evaluating them at non-interacting Kohn-Sham 
Green functions, one might expect a less stable functional 
to give rise to inferior response functions. And this is 
something which has to be thoroughly investigated. But 
it is clear that the Klein functional is much easier to 
evaluate and manipulate as compared to, e.g., the LW 
functional. 

All the $ functionals lead to a Dyson equation which 
has to be solved self-consistently for G. This is, in gen- 
eral, a very demanding task because of the complicated 
satellite structure inherent to any interacting Green func- 
tion. This severe complication is, however, circumvented 
by switching to TDDFT. 

Our approximations within TDDFT are just special 
cases of the variational functionals in which we re- 
strict the variational domain of the Green function to 
be all Green functions obtainable from a one-electron 
Schrodinger equation with a local multiplicative poten- 
tial - or vector potential in the case of current-DFT. 

We remark that this restriction on the variational 
freedom renders all the variational functionals density 
functionals iSSiS^ Given a density there is a local poten- 
tial which in a non-interacting system produces that den- 
sity. This potential produces the non-interacting Green 
function which we use to evaluate our functionals. Thus, 
the variational approach naturally generates different ap- 
proximations within DFT for static problems and within 
TDDFT for time-dependent problems or for the response 
functions of stationary problems. As we shall see, the 
exchange-correlation quantities depend on the choice of 
the action functional so that to every approximate Baym 
functional <f> correspond different approximate exchange- 
correlation potentials and kernels. 

Below, we discuss TDDFT and TD current-DFT (TD- 
CDFT) approximations in the framework of the Klein 
functional and of the LW functional. We also generalize 
the theory to ^' functionals and give some examples of 
approximations which we believe to be quite feasible to 
apply to realistic systems taking due account of the full 
electronic structure of one-body origin. 



III. TDDFT FROM THE KLEIN FUNCTIONAL 

Let Gg be the Green function of a non-interacting sys- 
tem of electrons exposed to the external, possibly time- 
dependent, potential V{rt). The Klein functional evalu- 
ated at Gs can then be regarded as a functional of V: 

iYkIV] = $[G,]-Tr {G.Gjj' - I + H-G;^)} ~tUn[Gs]. 

We could now directly use the stationary property of 
the Klein functional with respect to variations in the un- 
known one-body potential V in order to obtain an equa- 
tion for that potential. Because of the simplicity of a non- 
interacting Green function, however, the functional Ik 
can first be manipulated to acquire a physically appeal- 
ing form. This can, most easily, be seen in the static case 
elaborated below. The following equations are still valid 
in the case of time dependent problems and/or problems 
at elevated temperatures. This, however, requires some 
reinterpretations of standard DFT quantities like, e.g., 
Tg or C/h I which then become functionals on the Keldysh 
contour. For non-interacting Green functions the loga- 
rithm of the inverse of Gs is just the sum of the occupied 
eigenvalues contained in G^.^^ And the trace of GsG^^ — 1 
is just the integral of the particle density multiplied by 
the potentials V — w — Vb_- Expressing the eigenvalues of 
the one-electron Hamiltonian — V^/2 -t- y as expectation 
values then leads to. 



-mGs 



(6) 



Here, the quantity Ts[n\ is the well known functional for 
the kinetic energy of non-interacting electrons in the po- 
tential V - which produces the density n. Comparing 
now with standard DFT we see that the ^-functional pre- 
cisely plays the role of the exchange-correlation energy. 
This means that we may reuse standard DFT results and 
realize that the Klein functional is stationary when 



V 



5^ 

on 



(7) 



where the last term is the exchange-correlation poten- 
tial Wxc- Using the chain rule for differentiation we can 
rewrite Wxc as 



.(5$ 



(5$ 5Gs 5V 

JgIIvI^ 



(8) 



Remembering that the derivative of <& with respect to 
the Green function is just the self-energy E of our $ 
derivable theory and that the last factor is the inverse 
of the density response function Xs of non-interacting 
electrons, we finally arrive at the equation 



,(2, 3)A(3, 2; 1) d(23) = j Xs(l, 2)vW^) d2, (9) 
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where we have defined a generalized non-interacting re- 
sponse function A according to 



zA(2,3;l) 



SV{1) 



G42,1)G,(1,3), 



in terms of which we have Xs(l, 2) = A(l, 1; 2). 

Eq. (jnj is exactly the "linearized" form of the Shani- 
Schliiter (SS) equation, i.e., it can be obtained from the 
SS equationi'^* by replacing the interacting G with Gg, 
and S with = S[Gs]. Thus, the "Hnearized" SS equa- 
tion follows from a variational principle. We realize that 
the variational approach can be used to obtain succes- 
sively better approximations to the exchange-correlation 
potential v^c by making successively better approxima- 
tions to the functional (f>. In addition, the $- derivability 
and the variational property renders any approximation 
fully conserving. 

The full density response function expressed in the 
manner of TDDFT is^ 



X = Xs + Xsiv + Uc)x, /xc(l,2) 



dn{2) • 



The kernel /^c can now be obtained from one further 
variation with respect to the total potential V. The vari- 
ation of Vxc with respect to V can be expressed in terms 
of the exchange-correlation kernel /xc as 



6V{2) 



Ac(l,3)x«(3,2)d3. 



Thus, by varying Eq. we obtain the following equa- 
tion for /xc 



A. The exchange-only approximation 

Let us consider, for instance, the simplest approxima- 
tion for <f>, namely the Hartree-Fock approximation: 



from which we obtain 



Sx(l,2) = 

In this case 

'5Sx(l,2) 
SV{3) 



<5G(2,i; 



= w(l,2)G,(l,2). 



7;i;(l,2)G,(l,3)G,(3,2), 



and Eq. (|1UI) leads to the diagrammatic expansion in Fig. 
^ This approximation is also known as the exchange- only 
(EXO) approximation or sometimes the exact exchange 
(EXX) approximation. It has been evaluated earlier by 
several people4i^*SiL2i2ii£ 



^^^^ =-iv(l,2) 

FIG. 1: Exchange-correlation kernel in the EXX approxima- 
tion 



where 



X«(l,3)/xc(3,4)x.(4,2) d(34) 

A(l,3;2)A(3,4)G,(4,l)d(34) 
G,(1,3)A(3,4)A(4,1;2) d(34). 



A(l,2) = S,(l,2)-5(l,2Kc(l). 



(10) 



When the potential v^c has been obtained from Eq. 
the right-hand side of Eq. is a calculable expression 
for any given approximate <I> and no self-consistency is 
required. As an additional bonus, all occurring Green 
functions are non-interacting as opposed to interacting 
as one would have in most iterative schemes based on 
MBPT. (Consider, e.g., the response function of the 
time-dependent Hartree-Fock approximation.) 



B. The GW approximation within TDDFT 

Let us now go one level higher in the expansion 
of many-body perturbation theory and include all the 
screening diagrams. This is called the GW level with <I> 
given by 

<f[Gs] = ^Tr {Hl + tvGsGs}}. 

The expression of the self-energy in the GW approxima- 
tion becomes 

2,(1,2) =^G3(1,2)VF(1, 2); W=[l~vXsr'v. 

To calculate the variational derivative of the self-energy 
with respect to V we need to evaluate the change in the 
screened potential W. This can easily be constructed 
by observing that W~^ = v^^ — Xs and that SW/5V = 
~W[5W-'^/5V]W. The final result is displayed in Fig. H 
in terms of Feynman diagrams. 
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FIG. 2: Exchange-correlation kernel in the GW approxima- 
tion 



All Green functions are Kohn-Sham Green functions 
and all interactions are RPA screened interactions. This 
response function for which the "time-dependent GW 
(TDGW) response" would be a descriptive name, is 
presently too difficult to compute in real systems. Gel- 
dart and Taylor used it to investigate the effects of the 
static screening properties on the electron gasi^ It was 
used by Langreth and PerdewSS in the static long wave- 
length limit in order to extract gradient approximations 
for DFT. Richardson and AshcroftSSi have published an 
approximation to the TDGW response of the electron gas 
but only at imaginary frequencies. Another application 
of the TDGW^ response is due to Langreth et al.^^ and 
deals with Van der Waals forces. The TDGW response 
is generally believed to be very accurate but the compu- 
tation of the screened interaction is known to be a bottle 
neck in GW calculations on real solids. Unfortunately, 
the TDGW response contains two such complicated fac- 
tors (screened interactions). 



As in the case of only density variations, the simplic- 
ity of a non-interacting Green function again allows the 
Klein functional to be written in a much more convenient 
form. Using similar manipulations as in the beginning of 
Sec. mil we arrive at the expression 



(11) 



where we have used the normal convention to sum over 
repeated indices. Here, the functional for the non- 
interacting kinetic energy also depends on the physical 
current density j and not only on the density n. As 
before, the ^-functional plays the role of the exchange- 
correlation energy. We then realize that the functional 
Yk is stationary when 



V = w + V-ii + Vxc where 

^ ^cxt 



(12) 

on 

(5$ , , 

-w '''' 



Let us now focus on those system with a vanishing ex- 
ternal vector potential. Following the same steps as 
led to Eq. ©, i.e., the chain rule for differentiation, 
we obtain the "linearized" Sham-Schliiter equation of 
TDCDFTi^lS^ 



S,(2, 3)A^(3, 2; 1) d(23) = / Xs,p.(l, 2)A,c,.(2) d2 



(14) 

(Notice that A-^cf^ = (^xc, ^xc) in normal four- vector no- 
tation.) The generalized response function A^ appearing 
above is defined according to 

zAo(2, 3; 1) ^ = G.(2, 1)G.(1,3), (15) 



C. TDCDFT from the Klein Functional 

In TDCDFT the density n and the physical current 
density j are uniquely fixed by the external vector po- 
tential Acxt and the scalar potential wm^^ The coupling 
to the external fields is given by 



J 



cPr [Aext(rt) ■ ip(rt) + w{rt)n{rt% 



where w — w + A^^j/2 and jp is the paramagnetic current 
operator. According to our prescription we render the 
Klein functional a functional of = (n,j) by restricting 
the variational freedom of the Green functions to be all 
those Gs 's which are non- interacting and given by a local 
scalar potential and a vector-potential, — (V, A). It 
is convenient to consider the four- vector = (y. A), 
where V — V + as the independent variables since 

the four- vector density jp^^ = (n,jp) is coupled linearly 
to iu. 



zA(2, 3; 1) ^ = ^G.(2, l)[Vi - Vi]G.(l, 3), 

(16) 

and from this response function we obtain the Kohn- 
Sham density-density, current-density, and current- 
current response functions from the relation below 



Xs,p.(l,2) = 



M,(2)- 

The many-body response function 



(17) 



Xa<.(1,2) 



<^^oxt,iy(2) 



can be expressed in terms of the Kohn-Sham response 
function Xs.i-iu through the exchange-correlation kernel 

Xm^ ~ /Cs,^u ~1~ Xs,/jp(/xc,pcr ^ ^po ^ f^ocr)Xcri/ (1^) 
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where /xc,/i!y = SAxc,ii/Sj,j. In our variational scheme the 
equation for /xc is obtained from one further variation 
of Eq. 1)14(1 with respect to the Kohn-Sham potential 
Afj^. The corresponding response function X/if obeys 
the /-sum rule and Ward identities'^- since under a 
gauge transformation the scalar potential V and vector 
potential A change as in the exact CDFT, namely 
V ^ V + df /dt and A — > A + V/. In order to prove 
this property we change the external fields according to 
w — > w + df /dt, Acxt — > Acxt + V/ and we ask the 
question how the scalar potential V and vector potential 
A change at the stationary point. From Eqs. H12I13|I . 
it is straightforward to realize that V V + df/dt 
and A ^ A + V/ provided the exchange-correlation 
potentials change according to v^c Vy^c + A^c ■ V/ and 
Axe Axe- Taking into account that under this gauge 
transformation Gs(l,2) e~'f'^^^Gsil,2)e'f'-^\ it is a 
matter of very simple algebra to show that the linearized 



SS equation (|14|l is gauge invariant for any ^-derivable 
self-energy. 



D. The EXO within TDCDFT 

Let us consider, for instance, the exchange-only ap- 
proximation for the homogeneous electron gas. Extract- 
ing the time-ordered component of Eq. (|18|l and taking 
advantage of the translational invariance of the homoge- 
neous electron gas, we find"^"* 

^)Xs.aiy{q, oj) = V^,y{q, Lu) + Sf_,^{q, us) 

where all quantities are time-ordered and where V^j^ and 
Sp,i,, at zero temperature, are given by 



V^^(q,w)= J (fpd^kp^v{\p-k\)K 



^p+q/20p~q/2 



yp+q/20p-q/2 



<^ - £p+q/2 + £p-q/2 - «?7 ^ ~ '^p+q/2 + ^p ql2 + «77 ^ 

^k+q/2^k-q/2 ^k+q/2^k-q/2 
. ^ ~ ^k+q/2 + ^k-q/2 " *^ ^ " '^k+q/2 + ^k-q/2 + , 



(19) 



yp+q/20p-q/2 



^p+q/20p-q/2 



- £p+q/2 + + i'/)^ - £_p+q/2 + ^p-q/2 - ivY 

X {Sx(p + g/2)-Sx(p-g/2)}. 



(20) 



Here, we have denoted by p^, the four-dimensional 
vectors of components (l,p), (1,A;), while the Heaviside 
step functions 

9q — 6{eF — eq) and Oq = 1 — 6q 

contain the Fermi energy sf. 

In the large uj limit the sum Voo + 5*00 goes like 1 fuj^ 
and therefore Xoo = Xs.oo + 0{l/uj*). Since the residue 
of the second-order pole in Xs.oo only depends on the 
density, the approximated response function Xoo obeys 
the /-sum rule, as it should. 

E. Conservation laws 

As mentioned several times, the variational and 
<I>-derivable approach to TDDFT leads to density- 
functional approximations which preserve many physi- 
cal properties when the system is subject to external 
perturbations. Of course, TDDFT being a one-electron 
like theory with a multiplicative potential trivially obeys 
the continuity equation and thus particle conservation 



for any approximation to exchange and correlation. The 
conservation of other quantities will however depend on 
the choice of such approximations. 

In this subsection we will, as an example, show how 
momentum conservation follows from the general for- 
malism. In the one-electron like theory of TDDFT, the 
change of total momentum per unit time is simply given 
by J nV(w + Vh + i^xc)- The approximation to exchange 
and correlation is momentum conserving provided Wxc 
satisfies the zero force theorem. '^^ Designing exchange- 
correlation potentials that fulfill such a constraint is non 
trivial;^ and several well-known approximations are ac- 
tually not conserving. '^^•'^^ Below, we show that any ap- 
proximate Wxc generated by our variational approach is 
fully conserving. 

From Sec. Ill, we know that the change (5$ in the 
^-functional is just 

S^^i J wxc(l)<5n(l) dl (21) 

when we change the one-body potential from V to V+SV. 
In the variational approach a la Klein, Eq. (|21|) plays a 
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similar role as the Baym construction (5$ = Tr [E(5G]. In 
order to prove the conservation of the total momentum 
we have to shown that Wxc does not exert any force on 
the Kohn-Sham system. Let us shift all coordinates by 
the same time dependent infinitesimal vector d{t). The 
functional $ does not change since the interaction po- 
tential is invariant under translations. This implies that 



= (5$ = i J v^c{l)S{ti) ■ Vin(l) dl. (22) 

One partial integration and the fact that the vector 6{t) 
is arbitrary and independent of position gives 



J n{rt)'Vv^cirt) (fir = 0. 



(23) 



This means that there is no contribution from exchange 
and correlation to the total force applied to the system 
which is given by the classical expression F = — J n^w, 
as it should. 

The proof of momentum conservation in the presence 
of vector potentials and currents follows in a similar way 
from the corresponding result 



(24) 



which we obtained from the Klein functional. 



IV. LW FUNCTIONAL 

Let us now discuss the variational functional of Lut- 
tinger and Ward. From Eq. © we find 



iSYi 



LW 



Tr 



1 



We introduce the auxiliary Green function G according 
to 



G^ 



1 



Gh + GhSsG, 



i.e., G represents the first iteration toward the full self- 
consistent many-body Green function starting from the 
Kohn-Sham Green function Gg. Writing the total poten- 
tial V as 

y = w + Vh + Wxc 

and eliminating Gh between Gg and G, one obtains G = 
Gg + G[^g — Vxc]Gs, and thus 



LW 



SV{1) 



Tr G [Ss - Uxc] Gg 



5^g 



SV{1) 6V{1) 



In the Hartree-Fock approximation, E^. ~ Ex and Eq. 
(113 yields 



G(2, 3)Ex(3, 4)G,(4, 5)t;(2, 5)A(5, 2; 1) rf(2345) 
G(2, 3)i;xc(3)G,(3, 2)w(2, 4)xs(4, 1) d(234). 



This equation determines the exchange-correlation po- 
tential which, in turn, fixes the total potential V and 
then G and G. We also observe that this v^c is approxi- 
mately linear in the strength of the Coulomb interaction 
as is the ordinary exchange potential described in the 
context of the Klein functional. As in the case of the EXO 
approximation discussed above, the exchange correlation 
part of the response function /xc is obtained from one 
further variation with respect to the total potential V. 
The expression now becomes slightly more tedious but is 
still of the form O/xc = 5xc, where O is a known operator 
and 5xc is a calculable expression which does not involve 
/xc, see Fig. El 

f /xc 

CXX>O^CD0O= 





(25) 



FIG. 3: Exchange-correlation kernel calculated from the LW 
variational approach to TDDFT in the Hartree-Fock approx- 
imation. The thin lines represent the noninteracting Green 
function Gs while the tick lines represent the Green function 
G. 



The response function of the LW functional at the 
Hartree-Fock level is expected to be superior to that of 
the EXO due to the better stability of the LW functional 
as compared to the Klein functional. 

We also wish to emphasize that it made no difference 
to the resulting response function whether we used the 
LW functional or the Klein functional to derive it - pro- 
vided we allowed free variations of the Green function G. 
By restricting the possible choices of Green functions to 
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those produced by local potentials (TDDFT) the LW and 
Klein functionals give rise to different response functions 
at the same level of many-body perturbation theory. 



V. * FUNCTIONALS 

The main advantage of the ^E* functionals is that they 
give the possibility of using physical models for the 
screening, the calculation of which is actually a bottle- 
neck in practical applications. A word of caution is, how- 
ever, appropriate in this context. With model screened 
interactions (W's) there is usually no self-consistency 
with respect to W, a fact that might compromise the 
conserving property of the theory. The functionals 
have two independent arguments (G and W) resulting in 
terms linear in the deviation of the actual Green function 
from the self-consistent one when W is away from the 
value which renders the functional stationary. When the 
^' functionals are used to construct response functions of 
TDDFT the theory is, however, variational with respect 
to the one-body potential generating the non-interacting 
Green function - even when a model W is used. This fact 
actually restores several conserving properties although 
this has to be verified from case to case. For instance, 
choosing model W:s which, like the bare Coulomb in- 
teraction, are instantaneous and translationally invariant 
will clearly not spoil the conserving properties. 

The first 5* functional was constructed by ABL in 
lOOB.JJ* It has the appearance 



iYABL[G,W] = *[G,1^]-Tr {EG-ln(E-GH')} 
+ ^Tr {WP + \n{l-vP)}-iUH[G]. 

(26) 



The ABL functional is stationary with respect to varia- 
tions of G and W whenever G obeys Dyson's equation 
and W obeys the "contracted Bethe-Salpeter equation" , 
W — V + vPW. As for the <& functionals, the self-energy 
is obtained by taking the functional derivative of 4* with 
respect to G and the polarization P turns out to be the 
negative of twice the functional derivative of ^' with re- 
spect to W 



E(l,2) 



<5* 



(5G(2,1)' 



P(l,2) = -2 



(5* 



6W{2,iy 



Just as was the case for the pure $ functionals, we can 
add to any ^E" functional an arbitrary functional K[Q] of 
a quantity Q defined by Q = {y-^ ~ P[G, W]) ~ 1, 
with the properties 

i^[0] = ^[0] = 0. 

We then obtain a new ^' functional with the same sta- 
tionary point and the same value at the stationary point. 



An example of a simple functional obtained in this way 
is 

«>LWs[G, W\ - *[G, W^] - Tr {EG - In (E - G^) } 
+ iTr {Wv-^ - l + ln(W^w-i)} ^ iUY^\G\. 

Here, LWS stands for the simple version of the ^ func- 
tional based on the construction of Luttinger and Ward. 
As for the Klein version of <f>-derivable functionals, we 
expect that this functional is less stable than the orig- 
inal ABL functional of Eq. (|26|l . Also, it contains the 
same LW expression which led to the complicated result 
depicted in Fig. |21 and, for the moment, we deem this 
functional as less suitable for the construction of response 
functions. 

Another possibility is to, instead, add to Yabl the 
functional F\D\ of Eq. lO), thus obtaining the Klein ver- 
sion of ^f-derivable functionals 

*i"ABLK[G,VK] = ^\G,W\ 

- Tr {GGjj^ - 1 + In (-G^^) } 



-Tr {WP 



ln(l -t;P)} -iC/H[G]. 



Again, due to the simplicity of the "Klein" expression, 
we can here use the same manipulations as we applied 
to the original Klein functional in order to arrive at Eq. 
©. Thus, inserting the non-interacting Green function 
Gs into the functional Yablk, we then obtain 



wn 



(27) 



where 



-i^Gs^W] 

-Tr {WP[Gs, W] + In (1 - vP[Gs, W])] . 



Consequently, also this functional can be given the stan- 
dard DFT form and wc realize that it is stationary at the 
non-interacting Green function Gg produced by the local 
one-electron potential 



y = w + Vh + 



5E^ 

Sn 



(28) 



In fact, all functionals, be they of the $ or the 5" vari- 
ety, having the "Klein" form for their dependence on the 
external potential w have the nice property that the opti- 
mizing potential consists of the external potential w, the 
Hartree potential Vu , and the functional derivative of the 
exchange-correlation energy with respect to the density 
n. In Eq. the last derivative is calculated from the 

chain rule for differentiation giving the OPM-like equa- 
tion 



Xs(l,2)wxc(2) d2 - 



5E^c 
5V{1) 



A(3,2;l) 



E(2,3) + i j Aiy(4,5) 



(5P(5,4) 
(5Gs(3,2) 



d(45) 



d(23). 



The quantity AW is W - W = W - v/{l - vP) and we 
remind the reader that we are here allowed to use any 
model for W. In particular, we could choose W to be W, 
in which case our equation for the exchange-correlation 
potential Uxc reduces to the same expression as obtained 
from the "Klein version" of the $ formalism described in 
Sec. mil Furthermore, it is easily seen directly from its 
definition that the functional Yablk becomes indepen- 
dent of the choice of model W at the level of the RPA. 
Thus, at that level, this functional does not add any- 
thing to the previously discussed ^-derivable scheme at 
the same level (RPA). Being, for the moment, content 
with that level we will here not pursue the Iablk any 
further. 

Finally, by adding an appropriate choice for the func- 
tional K[Q], as discussed above, to the functional Yablk, 
we obtain the simplest functional Ikk of those discussed 
in the present work. We have 

iYKK = *-Tr {GGh^ - l + ln(-G-^)} 

+ {Wv-'^ -l + \n{Wv-^)] -iUn[G]. 

As before, we restrict the variational freedom of the 
Green functions to non- interacting ones (Gs) and differ- 
entiate Ikk with respect to the total, as it turns out, 
Kohn-Sham potential V producing the non-interacting 
Gs. Notice that W is an independent variable and does 
not depend on Gg or V . Only Gg depends on V . We 
obtain 

y"s,(2,3)A(3,2;l)d(23) = J Xsil,2)v^,{2) d2, (29) 

where v^c — V ~ w ^ Va, as before. We are now allowed 
to choose any appropriate but approximate W = Wo- 
Let us study the response function resulting from the 
functional Ykk at the GW level, meaning that 

* = ^Tr {G,G, W^o} . 

Consequently, in this approximation, we obtain the GWo 
self-energy 

= — - iGsWo. 
oGs 

One further variation of Eq. H29|l with respect to V gives 
an equation for /xc, whose diagrammatic representation 
is given in Fig. 01 Here, the Green functions are Kohn- 
Sham Green functions and the corresponding exchange- 
correlation potential is that which, to first order, repro- 
duces the density of a GW calculation with some model 
Wo- The screened interactions could be model interac- 
tions, e.g., Yukawa potentials, screened potentials within 
the RPA, or plasmon-pole approximations. The vertex 
diagram in Fig. ^should be calculated with Kohn-Sham 
Green functions. The sum of the first four self-energy di- 
agrams is essentially, to first order, equivalent to a single 
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AAAAA = -I W„ 




FIG. 4: Exchange-correlation kernel at the GW level in the 
"if scheme. 



polarization diagram calculated using GWo Green func- 
tions. 

If Wo is chosen to be a Yukawa potential, i.e., a stati- 
cally screened Coulomb interaction, this conserving GWo 
response function provides a justification of the work 
of Marini et. al^ provided the Green function with 
GW (RPA) shifted poles is close to the Green function 
of a Hartree-Fock calculation with a statically screened 
Coulomb potential. 

The use of a Yukawa potential for a screened inter- 
action also sheds some light on the work by Kim and 
Gorling^. They did exactly the TDDFT response func- 
tion of Fig. 2] in bulk silicon using the bare Coulomb 
interaction for W. They found that their calculated op- 
tical absorption spectrum was far from an experimental 
result unless they, somewhat artificially, cut down the 
range of the particle-hole interaction. 

It would actually be quite interesting to investigate 
the properties of the response function of Fig. 0] using 
a variety of screened interactions of the kind that pre- 
serves its conserving properties. An ongoing collabora- 
tive project^Si^ has precisely this objective in mind. 



VI. CONCLUSIONS AND OUTLOOKS 

In the present work we have proposed a new way of ob- 
taining approximations to current and density response 
functions of realistic systems. Our theory is based on 
the variational approach to many-body theory as previ- 
ously formulated by us and others. It gives the possi- 
bility to find successively better approximations to the 
effects of exchange and correlation in the framework 
of time-dependent density-functional or current density 
functional theory. The fact that the theory is formulated 
in the language of TDDFT makes it much easier, from a 
computational perspective, to apply to realistic systems 
as compared to standard MBPT. 

Improved approximations can be constructed in a sys- 
tematic way in the same sense as within MBPT. But 
physical intuition as to what physical processes are im- 
portant for any particular problem must be applied. The 
underlying variational approach to MBPT is not unique 
and many variational functionals can be constructed 
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leading to the same quality of approximation. Our 
method for improved approximations within TDDFT has 
the same feature. Different functionals have different 
variational accuracy meaning different sizes of the sec- 
ond order errors. In the present paper we have dis- 
cussed mainly two functionals - that due to Luttinger 
and Ward (LW) and that due to Klein (K). The former 
has proved to be more stable as compared to the lat- 
ter as far as concerns the calculation of total energies 
of a variety of systems ranging from those with very lo- 
calized electrons to those with itinerant electrons. This 
would suggest that the LW functional ought to be used 
also for the construction of response functions within 
TDDFT. In the present work, we have given the formulas 
for the exchange-correlation kernel of TDDFT resulting 
from both functionals. Sadly enough, we judge that of 
the supposedly better LW functional to be beyond our 
present computational facilities - even at a rather low 
level of approximation within MBPT. In order to demon- 
strate this point, we have given the diagrams representing 
the density response function resulting from the LW for- 
mulation within the exchange-only approximation. We 
would still like to draw the readers attention to the fact 
that the ambiguity in the choice of functionals can most 
likely be used to our advantage. But much more research 
is needed in order to see how this should be done. 

A very important feature of our variational approach 
to TDDFT is the fact it rehes on the $ or * dcriv- 
ability of the underlying approximation within MBPT. 
Combined with the variational property of the chosen 
functional, this leads to the preservation of many phys- 
ically important conservation laws and sum rules. And 
this is true regardless of the actual chosen level of ap- 
proximation within MBPT. This highly desirable feature 
is not guaranteed in other available approaches based 
on straight-forward diagrammatic expansions, iterative 
techniques, or decoupling schemes. For instance, in Ref. 

one develops a diagrammatic representation for the 
particular many-body perturbation scheme which starts 
from a zero-th order Hamiltonian which already gives 
the correct density. '^'^ Unfortunately, this technique suf- 
fers from the same basic drawbacks as ordinary MBPT - 
it is, in principle, divergent, summations must be carried 
to infinite order, and there is no guarantee for obtain- 
ing approximations which have certain desirable physi- 
cal properties automatically included. The same holds 
true for expansions which are based on iterating the so 
called Hedin equationg^^ using the screened interaction 
as the "small" parameter.'*^ As an example, we have, in 
the present work, demonstrated how the variational ap- 



proach leads to momentum conservation in the case of 
the Klein functional. 

It is worth while observing that the so called "lin- 
earized" Sham-Schliiter equation actually turns out to 
be a result of our variational approach starting from the 
Klein functional. But this is only true if the self-energy 
involved is a $- or ^P-derivable one. In that case, the 
resulting approximation for the response function is, of 
course, conserving. 

We also remark that the so called optimized potential 
method (0PM) and many generalizations thereof readily 
follows from the theory presented here. As an example, 
we have given the explicit formulas for the current density 
response of a homogeneous system within the exchange- 
only approximation. 

Even though we now have a systematic way of obtain- 
ing better response functions within TDDFT the expres- 
sions quickly become too complicated to be implemented 
in low-symmetry systems, especially when we want to in- 
clude all physically relevant processes. In this context, we 
advocate the use of the ^f-derivable theories which allows 
for the use of model screened interactions without loosing 
the important conserving properties. In this way, impor- 
tant physical effects like, e.g., a strong particle-hole inter- 
action can be incorporated without an excessive increase 
in the computational effort. One should, however, keep 
in mind that models for the screened interaction must 
possess certain symmetries related to the actual system 
in order for the conserving properties to be preserved. 

We have discussed the implementation of the theo- 
ries presented here with other research groups. One 
particularly promising approach is that which is based 
on the Klein functional and the ^'-formulation using a 
model screened interaction like, e.g., a statically screened 
Coulomb interaction as often used in the Bethe-Salpeter 
approach, or a simple plasmon-pole approximation. To- 
gether with our coUaborators^Si^ we hope to be able to 
present some numerical results within the near future. 
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